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Abstract 

We study cc color-octet contribution to B — * if)(vp')X. When this contribution 
is included, the theoretical predictions for the branching ratios become in 
much better agreement with the experiment. This mechanism also enhances 
the partial rate asymmetries by about a factor of five. The inclusive ip(ip') 
resulting from b — > d + gluon can have asymmetry around a few percent 
whereas those from b — ► s + gluon has it around 4 x 10 -4 . The asymmetry 
in the former modes should be observable, to a significance of 3a, with about 
(1 - 10) x 10 8 B mesons. 
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Recent experimental data from the Tevatron indicate that if only cc color-singlet contri- 
bution is included the ip' production rate at large transverse momentum predicted by QCD is 
about a factor of 30 below the experimental data. It has been shown by Braaten and Fleming 
|]J and Cho and Leibovich || that if the cc color-octet also contributes to the tp' production, 
the experimental data can be explained. Color-octet also has significant contribution to the 
ip production at the Tevatron ||. In this paper, we show that, if this mechanism is indeed 
the correct one, it also has important implications for decays of B mesons, especially 

for CP violating particle-antiparticle partial rate asymmetry. The point is that the penguin 
graph leads to an appreciable branching ratio (~ 10~ 2 ) for b — > s + gluon. Now, in a purely 
perturbative approach, the formation of ip{ip') from the gluon is severely suppressed. On the 
other hand, if Braaten et aVs mechanism tends to enhance the rate for the color-octet (cc 
or gluon) to form the ip(ip') then it can have important consequences for direct CP violation 
in inclusive B decays, via B — > ip{ip') + X. This mechanism also enhances the branching 
ratios in these decays so that they are much closer to the experimentally measured ones. 

In the SM the amplitudes for B decays are generated by the following effective Hamil- 
tonian: 

/-i 10 

H! ff = -^[ Vfb v; q ( Cl oi f + C2 ot f ) - E(K.K>r + v ch v^ct + v a vtf)of\ + h.c. , (i) 

V 2 i=3 

where the superscripts u, c, t indicate the internal quarks, / can be u or c quark, q can be 
d or s quark depending on if the decay is a AS = or AS = — 1 process. The operators Of 
are defined as 

0} x = qajpLfpfp-fLba , 0% = qjpLff-fLb , 

°h = q.ln Lb qtln L i R )<l' > = Qal»Lbpq , pj lx L(R)q' a , (2) 

3 3 
°7,9 = ifq'Ql^Lbq'YR{L)q' , 0l,io = ^q'Q a l^Lbpq'^^R{L)q' a , 

where R(L) = 1 + (— )7s, and q' is summed over u, d, s, and c. are the tree level 

and QCD corrected operators. O^s are the strong gluon induced penguin operators, and 
operators O 7 _i are due to 7 and Z exchange, and "box" diagrams at loop level. The WC's 
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c{ are defined at the scale of /x ~ m&. Although the WC's have been evaluated to the 
next-to-leading order in QCD H|J, we will use the leading order WC's to be consistent with 
the matrix elements which were evaluated using NRQCD to the leading order. We give the 
coefficients below for m t = 176 GeV, A4 = 0.2 GeV, and \x = m b = 5 GeV, 

ci = -0.249 , c 2 = 1.108 , c* = 0.0116 , c\ = -0.0249 , c\ = 0.0073 , c\ = -0.0300 , 

4 = 0.0011 , 4 = °- 0004 > 4 = -0.0092 , 4 = 0.0021 , 

cl$ = ~c u 4/N = P C JN , 4; 9 c = Pr, 4,^ = ; (3) 

where N is the number of colors. The leading contributions to P l se are given by: P* = 
(a s /8vr)c 2 (10/9 + G{m h n,q 2 )) and P\ = (a em /9vr)(iVci + c 2 )(10/9 + G{m h ^q 2 )). The 
function G(m, //, q 2 ) is given by 



G(m,iJ,,q ) = 4 / x(l - x)dxln ^ — , (4) 

Jo n z 



where q is the momentum carried by the virtual gluon in the penguin diagrams. When 
q 2 > 4m 2 , G(m, /1, q 2 ) becomes complex. In our calculation, we will use m u = 5 MeV and 
m c = 1.35 GeV. Although there is considerable uncertainty in the effective values of these 
masses, especially so for m u , our results are rather insensitive to the numerical values of 
these masses. 

Using factorization approximation, we have 

M(B - J/ijX s ) = %A X < Xl\s lfl (l - l5 )b\B >< 1>W)\cf(\ - ^clX 1 > 
+ 2A 8 < X 8 s \s^(l - j 5 )T a b\B >< - 7 5 )T a c|A 8 > 

^tu ^tu ^tu ^tu 

A, = V cb V; s (c 2 - c? - cg» - ct - cTo) - V tb V*(cf + cf + of + c%) . (5) 

where cf = c? — c" and 4" = 4~4' + = X^ + X s = X s . The term proportional to Ai 
is the color-singlet amplitude and the term proportional to A 8 is the color-octet amplitude. 



If the color-octet contribution is neglected, the branching ratios for B — > ip(ip')X are 
too small compared with the experimental values. In order to reproduce the experimental 
data, the number of colors N is traditionally treated as a free parameter to parametrize 
the non-factorizable and the color-octet contributions. The effective number of colors N is 
then determined from B — > %p(ip')X to be close to 2 [[J. This does not really identify where 
the new contributions come from. However, if the color-octet effects identified above have 
significant contributions, one may not need to treat N as a free parameter; N = 3 as given by 
QCD may work fine. The color octet mechanism with N = 3 indeed improves the situation 
significantly. This has been pointed out by Ko, Lee and Song ||. Our calculations confirm 
their results. In their work the penguin contributions are not included. When penguin 
contributions are included, they have important implications for direct CP violation in 
particle-antiparticle partial rate asymmetries in these decays although their effect on the 
absolute rates is minimal. 

In order to obtain non-zero partial rate asymmetry, it is necessary to have non-zero CP 
violating phases and CP conserving strong phases due to the final state rescattering for 
different amplitudes. In the above case, the CP violating phases are provided by the phases 
in the CKM matrix elements V c t,V* s and V t bV* s . We will use the Wolfenstein parametrization 
such that the element V u b is given by |Vub|e _i7 . For the strong phases, we will use the 
phases generated at the quark level by appealing to quark-hadron duality. This should 
be at least a good indication for the size of the phases j7|. Naively the strong phases are 
generated by exchanging u and c quarks in the loop. However, CPT theorem dictates that 
the phases generated by the c quark in the loop not to contribute to the rate asymmetry for 

. Both the strong penguin WC's c^ i56 and the electroweak 
penguin WC's c" 8 9 10 contribute to the strong phases. It is, however, interesting to note that 
if the color-octet amplitude is neglected, the strong penguin do not generate non-zero strong 
phases because the combination c% +c^/N + c^ + c$/N in A 1 is identically zero as can be seen 
from Eq|| The leading non-zero strong phases are then generated by electroweak penguin 
and are therefore small. If these phases are the only ones, the partial rate asymmetry is 
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the production of ip{ip') [B 



predicted to be very small as shown in Figure 1 flO] . Here we have used N = 2 since no 
color-octet contributions have so far been included as discussed before. If the color-octet 
contribution turns out to be significant, the situation can become dramatically different. 
Including the color-octet contribution, we have 

\M(B - </#')X s )| 2 = G 2 F Tr(f s + ro.M/ft + m b )Y(l - 7s) 



x (-fiv + 



x [i^r+4iv <Q8 !lm (3gl)> ], (6) 



where the operators Of ^ \ 3 Si) are defined in Ref. pjjjll 

From eq|], the branching ratios and the CP violating partial rate asymmetries can be 
easily calculated. To quantitatively assess the importance of the color-octet contribution 
we recall that it has been shown by Ref. that even with a small matrix element for 
color-octet to produce a ip{ip')-, the experimental data from the Tevatron can be understood. 
Fitting the experimental data from the Tevatron, Cho and Leibovich obtain [0 

< OtCSi) > = 1.2 x KT 2 GeV 3 , 

< Of ( 3 Si) > = 7.3 x KT 3 GeV 3 . (7) 
The color-singlet matrix elements determined from leptonic decays of ip and ip' are 

< Of ( 3 5i) >= 1.32GeV 3 , 

< Of' ( 3 Si) >=0.53GeV 3 . (8) 

Without the color-octet contributions, and with N = 3, the branching ratios predicted, 
are several times smaller than the experimental values: Br(B ipX) = (0.8 ± 0.08)%, and 
Br(B —>■ ip'X) = (0.34 ± 0.05)%. The inclusion of the color-octet contributions improves 
the situation with the branching ratios now predicted to be: Br(B — > ipX) = 0.54%, and 
Br(B — > ip'X) = 0.25%, with N = 3. These numbers are in good agreement with Ko et al 
0] and they are also much closer to the experimental values. 



Unfortunately the results are very sensitive to c\ and c%- The dominant color-singlet con- 
tributions are from operators 0\ 2 which are proportional to Ci+c^/N. There is an accidental 
cancellation here. Had one used the next-to-leading coefficient for ci,2, the cancellation is 
even more severe. If one adjusts the scale \i and A4 for the leading coefficient, one can get 
larger values for the branching ratios. There are other uncertainties in the evaluation of the 
branching ratios, namely, there are more operators which may contribute to the branching 
ratio. For example, there may be contributions from Ofg ^(^o)- The value for < 0\QS§) > 
is expected to be much smaller than < 0\( 3 So) >■ Its contribution to the branching ratios 
is expected to be small. The contributions from Os^Sq) may be potentially large because 
< 8 ( 1 S'o) > may be not too much smaller than < O 8 ( 3 S , ) >. If this is indeed the case, the 
experimental branching ratios can be easily reproduced. 

When the color-octet contributions are included the strong penguin also generate strong 
phases through the coefficient A$ in Eq. |5|. These phases are much larger than the ones 
generated by the electroweak penguins, and therefore much larger partial rate asymmetries 
result. The results are shown in Figure 2. In the figure, we used q 2 = vr?^,^ because the 
ip{ip') carries most of the momentum from the virtual gluon. We also set a s at q 2 = m^w) 
and the corresponding value a,(m^u,)) = 0.27 for A 4 = 0.2 GeV. If larger a s is used, the 
asymmetries become bigger. The asymmetry for B — > ip'X s is slightly larger than that 
for B — > ipX s . This is due to the fact that the ratio of the color-octet matrix element to 
color-singlet is larger for the if) r . We also considered the contribution from dipole penguin 
operators, O n = {g s /l^ 2 )m b sa^RT a bG^ and 12 = {e/lQii 2 )m b sa^RbF% u . Here G>" 
and F^ v are the gluon and photon field strengths, respectively. It has been shown that 
the operator On can enhance certain B decay branching ratios by as much as 30% ||14|| . 
However, its contribution to B — > ip(ip')X s branching ratio is less than 10 -4 and to partial 
rate asymmetry is less than 10 -5 . O12 contributions are even smaller. We remark that 
even if the operator Os( 1 5'o) contributes significantly to the branching ratios, it will not 
introduce new strong phases in the amplitude because its contributions are proportional to 
(c" — Cg ) which generate vanishing absorptive amplitudes. And therefore, it will not 
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affect the asymmetries evaluated here. 

It is clear from comparison of Fig.l and Fig. 2 that inclusion of the color-octet enhances 
the asymmetries. The asymmetries are bigger by about a factor of five. At present the CP 
violating phase 7 is not well determined; sin 7 can vary from 0.1 to 1. If we use the best fit 
value from the experimental data, 7 is about 70° [I5|. With this value, the rate asymmetries 



for B — > ip{ip') are about 4 x 10~ 4 (6 x 10~ 4 ). In order to observe the asymmetries at the 
3cr level in B — > ipX s and B — > ip'X s , one would need about 4 x 10 9 B decays. This number 
does not include any factor(s) for experimental efficiencies. So the number of B's needed is 
likely to be even higher depending on the specific final states of the ip(ip') that are accessible. 

The situation with B — > ■?/> (?/>') A^ (A^ denotes states without strangeness number) is 
better. The analysis is similar to B — > ip(if)')X s case. One only needs to change the 
relevant CKM matrix elements V c t,V* s and V t bV t * to V c bV* d and V t bV t * d , respectively in Eq. ||. 
The results are shown in Figure 3. The asymmetries can be as large as 1.5%. They are 
about 1% (1.5%) for B — > ip(ip')X s with 7 = 70°. It is interesting to observe that these 
asymmetries are similar to those obtained in Ref. ]ll]] by considering absorptive contribution 
from rescattering of cc color-octet states. For the X d final state the branching ratios are 
smaller by a factor of |V^/V^ S | 2 compared with B — > ip(ip')X s . Therefore to observe the 
asymmetries in B — > if)(if)')Xd at the 3a level, one would need about 1 x 10 8 B decays. 
Assuming an effective efficiency of 0.2 (i.e. including the branching ratio into some specific 
final state(s)) would make the actual number be more like 5 x 10 8 . 

The number of B's needed is clearly rather large so that even B factories may have a 
difficult time. On the other hand the ip and ip' tend to give distinctive signal which may 
even be accessible in a hadronic environment with a B-detector. 

In the above, we have considered the parton level processes, b — > s(d)ifj(tfj'). If one 
considers the parton level processes, b — > s(d)Y with Y being ip and other cc states which 
can materialize into tp, and do not isolate each individual Y and let it decay into ip, then the 
hadron level process, B — > ipX will have more sources to provide strong phases. In addition 
to the parton level strong phases discussed above, there is also the possibility of generating 



strong phases from resonant effects [[H]. The partial rate asymmetries may be even larger 
than what we obtained here. We will discuss the results from these mechanisms elsewhere. 

One might think that the same mechanism will also enhance the rate asymmetry in 
B — > r] c X s (Xd). It turns out that this is not true here if only Oi i s( 3 5'o) operators are included. 
In this case, for the same reason as for B — > ip(if)')X, the color-singlet only generates strong 
phases through electroweak penguins. However, the inclusion of color-octet contributions 
will not improve the situation because the analogous color-octet parameter A 8 (r] c ) is different 
than the Ag parameter in Eq.[5| The strong penguin contribution in Ag(i] c ) is proportional to 
^ c «(*) c _ C gW c ) wri ich generates vanishing strong phase too. However, if the operator Og^So) 
contributes significantly, strong phases will be generated by the strong penguins because the 
contributions are then proportional to (c^ c + Cq^ c ). There may be sizable rate asymmetries 
which also need further study. 
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FIGURES 




30 60 90 120 150 180 
FIG. 1. The partial rate asymmetry for B — > ip(i/j f )X s without color-octet contribution with 
\V c b\ = 0.04, iKtb/Fcfel = 0.08 and \V US \ = 0.22. The vertical axis is the asymmetry and the 
horizontal axis is the value in degree for the phase angle 7 in the Wolfenstein parametrization. 




FIG. 2. The partial rate asymmetry for B — > i/j(ip')X s . The solid and dashed lines are for 
B — > ipX s , and B — > ij)'X s , respectively. 
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